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A method of monitoring a magnetic field drift of a magnetic resonance imaging apparatus 



FIELD OF THE INVENTION 

The present invention is related to the field of magnetic resonance imaging 
(MRI), and more particularly without limitation to magnetic field drift compensation. 

5 BACKGROUND AND PRIOR ART 

Prior art magnetic resonance imaging apparatus 9 generate a uniform magnetic 
field which is also referred to as the main field or the polarising field Bo. The purpose of the 
magnetic field is that the individual magnetic moments of the spins in the tissue to be 
visualised attempt to align with this polarising field, but precess about it in random order at a 
10 characteristic Larmor frequency which is determined by the gyromagnetic constant of the 
spins and the polarising magnetic field B 0 . 

However during operation of a magnetic resonance apparatus the magnetic 
field can vary over time due to various reasons. For example when the magnetic resonance 
apparatus is under heavy duly the passive iron shim pieces of the coils which generate the 
1 5 magnetic field heat up which leads to a fluctuation of the magnetic field 

U.S. Patent No. 6,294,9 13bl shows a method for compensation of variations 
in the polarising magnetic field during magnetic resonance imaging. Monitor signals are 
acquired in an interleaved maimer during a scan with the MRI system. Frequency changes 
caused by variations in the polarising magnetic field Bo are measured using the monitor 
20 signals, and these measured frequency changes are employed to compensate image data 

acquired during the scan. This compensation is achieved by changing the frequency of the RF 
transmitter and receiver to offset the effects of changes in B 0 . A disadvantage of this 
compensation method is that special monitor signals are required for the compensation. 

It is therefore an object of the present invention to provide for an improved 
25 method of monitoring and compensating a magnetic field drift 

SUMMARY OF THE INVENTION 

The present invention provides for a method of monitoring a magnetic field 
drift of a magnetic resonance imaging apparatus which does not require special monitor 
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signals or monitor sensors. Rather the method of the invention can be performed on the basis 
of the magnetic resonance signals which are acquired during magnetic resonance image data 
acquisition, such as by single shot EPI or by a gradient echo sequence. 

In essence the phases of at least two magnetic resonance signals are acquired 
5 an echo time after the corresponding RF excitations. This corresponds to the central k-space 
line which has frequency encoding but no phase encoding. This central k-space is usually 
designated as k y = 0. The difference of two consecutive phase measurements which are 
acquired at a certain time interval provides the shift of the resonance frequency. This enables 
monitoring of the shift of the resonance frequency and compensation of the magnetic field 
10 drift. 

In accordance with a further preferred embodiment of the invention a single 

shot EPI method is used f or the ma gnetic resonanc e data acquisition. Single shot EP I has the 

advantage that every data acquisition contains the k y = 0 line such that for each data 
acquisition the required phase information can be obtained. 
15 In accordance with a further preferred embodiment of the invention a gradient 

echo sequence is used for the magnetic resonance image data acquisition. During such a 
gradient echo sequence the k-space is scanned along an arbitrary trajectory which typically 
also contains the k-space line k y = 0. For example a full gradient echo sequence has 256 data 

agquMtionS-Qnej)f which is re presentative of k y = 0, The data ac quisition for , k y = 0 enables 

20 to obtain the required phase information once for each complete gradient echo sequence. 

In accordance with a further preferred embodiment of the invention extra 
scans along k y = 0 are performed during a gradient echo sequence in order to obtain multiple 
phases during the complete sequence. This enables to monitor a magnetic field drift which 
occurs during a single gradient echo sequence. 

25 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the following preferred embodiments of the invention will be described in 
greater detail by making reference to the drawings in which: 

Figure 1 shows a block diagram of a magnetic resonance apparatus with 
30 magnetic field drift compensation, 

Figure 2 is illustrative of the phase determination for a single shot EPI, 

Figure 3 is illustrative of a flow chart for monitoring and compensating the 
magnetic field drift in the case of single shot EPI, and 



PHNL03 1 073EPP 

3 16.09.2003 
Figure 4 is illustrative of a flow chart for magnetic field monitoring and 
compensation if a gradient echo sequence is employed 

DETAILED DESCRIPTION 
5 Figure 1 shows a magnetic resonance device 1 which includes a first magnet 

system 2 for generating a steady magnetic field, and also several gradient coils 3 for 
generating additional magnetic fields having a gradient in the X, Y, Z directions. The Z 
direction of the co-ordinate system shown corresponds to the direction of the steady magnetic 
field in the magnet system 2 by convention. The measuring co-ordinate system x, y, z to be 

10 used can be chosen independently of the X, Y, Z system shown in FIG. 3. The gradient coils 
are fed by a power supply unit 4. An RF transmitter coil 5 serves to generate RF magnetic 
excitation pulses and is connected to an RF transmitter and modulator 6. 

A receiver coil is used to receive the magnetic resonance signal generated by 
the RF field in the object 7 to be examined, for example a human or animal body. This coil 

15 may be the same coil as the RF transmitter coil 5. Furthermore, the magnet system 2 encloses 
an examination space which is large enough to accommodate a part of the body 7 to be 
examined. The RF coil 5 is arranged around or on the part of the body 7 to be examined in 
this examination space. The RF transmitter coil 5 is connected to a signal amplifier and 
demodulation unit 10 via a transmission/reception circuit 9. 

20 The control unit 1 1 controls the RF transmitter and modulator 6 and the power 

supply unit 4 so as to generate special pulse sequences which contain RF pulses and 
gradients. The phase and amplitude obtained from the demodulation unit 10 are applied to a 
processing unit 12. The processing unit 12 processes the presented signal values (also 
referred to as k-space) so as to form an image by transformation. This image can be 

25 visualized, for example by means of a monitor 1 3 . 

Processing unit 12 determines a shift Af of the resonance frequency on the 
basis of the phase information obtained from the demodulation unit and provides the shift 4f 
to the control unit 1 1 for compensation of the shift. 

The acquisition of the required phase information is explained in more detail 

30 in the following by making reference to figure 2. First an RF magnetic excitation pulse 200 
is generated by the transmitter coil (cf. transmitter coil 5 of figure 1). Further a slice selecting 
gradient Gz pulse 202, a bipolar switched gradient 204, and a phase encoding gradient 206 
are applied. The application of these gradients 202, 204 and 206 for single shot EPI magnetic 
resonance data acquisition is as such known from the prior art 
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The single shot EPI provides measurement signal 208, The peak of 
measurement signal 208 is the echo time TE after the RF magnetic excitation pulse 200, The 
phase of measurement signal 208 at time TE with respect to the RF magnetic excitation pulse 
200 is determined for the purpose of monitoring a drift of the magnetic field and a 
5 corresponding drift of the resonance frequency. In order to determine the resonance 

frequency drift, if any, the phase of the measurement signal 208 is determined again at a 
subsequent single shot EPI. The phase difference between the single shot EPIs enables to 
calculate the shift of the resonance frequency. This will be explained in greater detail by 
making reference to figure 3. 
10 The phase information at time TE can be obtained from the measurement 

signal 208 in the time domain or in the frequency domain which requires to perform a Fourier 

transfor mation on measureme nt si gnal 208. The requir ed phase infor mation corresponds to a 

central line of the k-space which is usually referred to as k y = 0, i.e. a scan with frequency 
encoding but without phase encoding. The feet that there is no phase encoding enables to use 
1 5 drifts of the phase information obtained from the measurement signal 208 for determination 
of the resonance frequency shift. 

Figure 3 shows a flow chart for monitoring and compensating of a magnetic 
field drift in the case of single shot EPI. In step 300 the index d is set to zero. The shift of the 

resonanceJfrequency A/1 is alsojsffifcj&jzm^ 

20 In step 302 the first single shot EPI is performed and the measurement of the 

phase <]> a (d) is performed as explained with respect to figure 2. 

The time interval of TE after the first single shot EPI is performed in step 302 
the second single shot EPI is performed in step 304. Again the phase <|> a is measured. 

In step 306 a shift of the resonance frequency is calculated by calculating 
25 the difference between the phases (]> a which have been acquired in steps 302 and 304 and by 
dividing the difference by 2% x the time interval between the consecutive single shot EPIs. In 
this case this time interval is the echo time TE. 

In step 308 the absolute value of the shift Af is compared to a threshold value. 
If the threshold value is surpassed a drift compensation is performed in step 310. Preferably 
30 this is done by adjusting the frequency of the RF magnetic excitation pulses to the shifted 
resonance frequency. 
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If the shift Af does not surpass the threshold no drift compensation is 
required. In step 3 12 the index d is incremented and the control goes back to step 304 for 
continuous monitoring the shift Af . 

Figure 4 shows an alternative embodiment which uses a gradient echo 
5 sequence rather than single shot EPI. Step 400 corresponds to step 300 of figure 3. In step 
402 the magnetic resonance data acquisition is started by starting to scan the k-space along a 
given trajectory. The trajectory may or may not cover the central line of the k space which is 
k y = 0. 

In step 404 the scan of the k-space along the trajectory is interrupted in order 
10 to perform a data acquisition for the central k-space line with frequency encoding but no 
phase encoding, i.e. k y = 0. This way the phase information (j) a is obtained 

In step 406 the scan along the k-space trajectory is continued. After a certain 
time interval the scan along the k-space trajectory is interrupted again in order to perform 
another data acquisition for k y = 0 to obtain another phase information <j> a . 
15 In step 410 the shift of the resonance frequency is calculated by calculating the 

difference between the phases determined in steps 404 and 408 and by dividing the difference 
by27C TE. 

In step 412 the absolute value of the frequency shift Af is compared to a 
threshold value. If the threshold value is surpassed drift compensation is performed in step 

20 414. In step 416 the index d is incremented and the control goes back to step 406 in order to 
continue the scan along the k-space trajectory. This procedure continues during the entire 
gradient echo sequence for continuously monitoring the shift of the resonance frequency and 
compensating the drift of the magnetic field if necessary. 

The RF pulses of this modified gradient echo MRI sequence are applied at a 

25 single frequency (for a single slice); the phase measured at the echo (without encoding 

gradients) reflects only the phase difference accumulated between the RF pulse and the TE in 
case that the NMR resonant frequency (due to B0) is not equal to the excitation frequency 
(excluding susceptibility and chemical shift effects). When the magnet B0 is such that the 
NMR resonant frequency is the same as the excitation frequency of the RF pulse then the 

30 phase error accumulated between RF pulse and TE will be zero or constant from TR to TR. 

When the B0 is such that it corresponds to a different resonant frequency than 
the one the RF pulse is exciting, then the phase error accumulated between RF pulse and TE 
will be a finite value (proportional to the difference and TE). At the peak of each RF pulse, 
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all magnetization is in phase, irrespective of the BO. Only after this time does the BO have an 
effect on the magnetization. 

So, each RF pulse behaves like a phase reset (as far as the fl) measurement 
process is concerned). The relative phase error (relative to zero at the RF pulse peak) at the 
TE increases as BO moves further from the RF excitation frequency. Other mechanisms 
(chemical shift, susceptibility) can contribute to the phase error measured at the TE. 
Assuming that these other contributors to phase error are constant from TR to TR, the change 
in BO can be determined by calculating the difference in phase errors between data acquired 
form two different TR's. 
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CLAIMS: 



1 . A method of monitoring a magnetic field drift of a magnetic resonance 

imaging apparatus, the method comprising the steps of: 

- performing a first data acquisition by a first magnetic resonance signal being caused 
by a first excitation, 

5 - determining a first phase of the first magnetic resonance signal an echo time after the 

first excitation, 

- performing a second data acquisition by a second magnetic resonance signal a time 
interval after the first data acquisition, the second magnetic resonance signal being 
caused by a second excitation, 

10 - determining a second phase of the second magnetic resonance signal the echo time 

after the second excitation, 

detennining a shift of a resonance frequency based on a difference of the first and 
second phases. 

15 2. The method of claim 1 , whereby the first and second data acquisition are 

performed using a signal shot EPI method. 

3. The method of claim 1, whereby the first and second data acquisitions are 
performed by means of a gradient echo sequence method. 

20 

4. The method of claim 3, whereby a k-space is scanned and second data 
acquisitions are performed intermittently to determine second phases in order to continuously 
monitor the shift of the resonance frequency. 

25 5. The method of claim 4, whereby the second data acquisitions are performed 

after fixed time intervals. 
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6. The method of any one of the preceding claims 1 to 5, further comprising 

compensating the magnetic field drift by changing the frequency of the excitation in 
accordance with the shift of the resonance frequency. 

5 7. The method of any one of the preceding claims 1 to 6, further comprising 

compensating the magnetic field drift by adjusting the magnetic field 

8. The method of any one of the preceding claims 1 to 7, further comprising 

comparing the shift of the resonance frequency to a threshold value and compensating the 
10 magnetic field drift if the threshold value is surpassed. 

_ The method of any one of t he preced ing claims 1 _to_8 y ^^^byj^jtotmd 

second phases are determined in the time domain. 

15 10. The method of any one of the preceding claims 1 to 8, further comprising 

performing a Fourier transformation of the first and second magnetic resonance signals and 
determining the first and second phases in the frequency domain. 

1L Aj*mffllfer_pm^ digital storage medium, fox 

20 monitoring a magnetic field drift of a magnetic resonance imaging apparatus, the computer 
program product comprising program means being adapted to perform the steps of: 

- detennining a first phase of a first magnetic resonance signal, an echo time after a 
first excitation, 

- determining a second phase of a second magnetic resonance signal the echo time after 
25 a second excitation, whereby the second magnetic resonance signal is acquired a time 

interval after the first magnetic resonance signal, 

calculating a shift of a resonance frequency based on a difference of the first and 
second phases. 



30 12. The computer program product of claim 1 1 , the program means being adapted 

to continuously monitor the shift of the resonance frequency. 
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1 3 . The computer program product of claim 1 1 or 1 2, the program means being 

adapted to control an excitation synthesiser in accordance with the shift of the resonance 
frequency. 

5 14. The computer program product of claim 1 1 , 1 2 or 1 3, the program means 

being adapted to control the magnetic field in accordance with the shift of the resonance 
frequency. 

15. A magnetic resonance imaging apparatus comprising processing means (12) 
10 for determining a first phase of a first magnetic resonance signal an echo time (TE) after a 

first excitation, for determining a second phase of a second magnetic resonance signal the 
echo time after a second excitation, the second magnetic resonance signal being acquired a 
time interval after the first magnetic resonance signal, and for calculating a shift of a 
resonance frequency based on a difference of the first and second phases and the time 
15 interval. 

1 6. The magnetic resonance imaging apparatus of claim 1 5 having display means 
(13) for displaying of the shift of the resonance frequency. 

20 17. The magnetic resonance imaging apparatus of claims 15 or 16 further 

comprising control means (1 1) for controlling the generation of the excitations in accordance 
with the shift of the resonance frequency. 

1 8. The magnetic resonance imaging apparatus of claims 15, 16 or 17, further 

25 comprising control means for controlling of the magnetic field in accordance with the shift of 
the resonance frequency. 
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ABSTRACT: 



The invention enables to monitor a magnetic field drift of a magnetic 
resonance imaging apparatus on the basis of the magnetic resonance signals, which are 
acquired during magnetic resonance image data acquisition, such as by single shot EPI or by 
a gradient echo sequence. The phases of at least two magnetic resonance signals are acquired 
5 an echo time after the corresponding RF excitations. This corresponds to the central k-space 
line, which has frequency encoding but no phase encoding. The difference of two 
consecutive phase measurements, which are acquired at a certain time interval provides the 
shift of the resonance frequency. This enables monitoring of the shift of the resonance 
frequency and compensation of the magnetic field drift. 

10 

Figure 3 
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